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Surface  impact  simulations  of  helium  nanodroplets 
Final  performance  report:  AFRL/RWK  grant  FA8651-11-1-0005 

Robert  J.  Hinde,  Department  of  Chemistry,  University  of  Tennessee 


1.  Introduction 

When  a  helium  (He)  nanodroplet  passes  through  a  pick-up  chamber  that  contains 
a  low-density  atomic  or  molecular  gas,  the  atoms  or  molecules  that  impinge  upon  the 
droplet’s  surface  frequently  stick  to  the  droplet  and  migrate  into  the  droplet’s  interior.  If 
several  of  these  dopant  species  dissolve  in  a  single  nanodroplet,  they  may  coalesce  inside 
the  droplet  and  form  a  polyatomic  aggregate.  The  low  temperature  (T  0.37  K)  inside  He 
droplets  can  stabilize  these  aggregates  in  relatively  high-energy  structures  that  represent 
local,  but  not  global,  minima  on  the  aggregate’s  potential  energy  surface.  Helium  droplets 
could  thus  potentially  hud  applications  as  nanoscale,  cryogenic  “reaction  chambers”  for 
the  preparation  of  exotic  molecular  species  that  have  unusual  chemical  properties.  Such 
an  approach  could  be  used  in  the  near  future  to  create,  in  an  atom-by-atom  fashion, 
nanostructured  materials  that  serve  as  next-generation  catalysts,  propellants,  or  munitions. 

For  synthetic  chemists  to  fully  realize  this  potential,  they  must  devise  methods  for 
removing  the  solvated  aggregates  from  the  He  droplet  in  a  nondestructive  manner.  One 
simple  approach  for  doing  this  might  be  to  gently  “land”  the  doped  He  droplet  on  a  solid 
surface;  the  hope  here  is  both  that  the  droplet’s  bulk  would  cushion  the  impact  of  the 
dopant  aggregate  with  the  underlying  surface  and  that  the  droplet  would  fragment  upon 
impact,  with  individual  He  atoms  or  small  He^  clusters  carrying  away  the  center-of-mass 
kinetic  energy  of  the  doped  droplet. 

We  summarize  here  progress  made  on  a  two-pronged  program  of  research,  supported 
by  AFRL,  that  uses  theoretical  and  computational  methods  to  investigate  both  (1)  the 
formation  of  small  metal  atom  aggregates  inside  He  nanodroplets  and  (2)  the  desolvation 
and  isolation  of  these  aggregates  through  gentle  surface  impacts,  or  “soft  landings” . 

The  hrst  prong  of  work  focused  on  the  development  of  a  highly  parallel  code  suite, 
QDROP,  that  simulates  the  structure  of  neat  and  atom-doped  He  droplets  using  variational 
path  integral  (VPI)  quantum  Monte  Carlo  methods.  We  have  used  QDROP  to  study  He 
droplets  with  as  many  as  N  =  1000  atoms;  we  believe  these  are  the  largest  droplets  studied 
to  date  using  the  VPI  quantum  Monte  Carlo  approach.  Our  results  provide  the  hrst 
detailed  atomistic  picture  of  the  dense  solvation  shells  that  surround  atomic  impurities 
solvated  in  He  droplets;  it  has  been  speculated  that  if  two  Mg  atoms  or  two  A1  atoms 
dissolve  in  a  single  He  nanodroplet,  these  dense  shells  may  impede  recomination  of  the  two 
metal  atoms. 
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The  second  prong  of  work  took  initial  steps  toward  developing  a  “dressed  classical 
dynamics”  code  for  simnlating  the  real-time  dynamics  of  neat  and  doped  He  droplets.  In 
this  approach,  individnal  He  atoms  are  treated  as  classical  particles  that  are  “dressed”  with 
a  single-particle  wavefnnction  that  represents  each  atom’s  delocalization  in  the  droplet. 
Two  He  atoms  in  the  simnlation  interact  via  a  potential  energy  function  that  is  obtained  by 
averaging  the  bare  He-He  interaction  over  the  two  atoms’  probability  density  distribntions. 

2.  Summary  of  work  completed 

2.1  Quantum  Monte  Carlo  studies  of  doped  helium  droplets 

Mg-doped  and  Al-doped  He  nanodroplets  have  recently  been  investigated  in  several 
experimental  laboratories  [1-3] .  These  nanodroplets  are  prepared  in  pick-np  experiments, 
in  which  large  He  droplets  pass  throngh  a  pick-np  cell  that  contains  a  low  pressnre  of  Mg  or 
A1  vapor;  by  adjnsting  the  metal  atom  vapor  pressnre,  the  typical  nnmber  of  metal  atoms 
picked  np  by  the  droplet  can  be  controlled.  The  experiments  snggest  that  when  a  single 
droplet  captnres  several  Mg  atoms  or  several  A1  atoms,  the  captnred  atoms  do  not  form  a 
strongly  bonnd  metal  clnster  inside  the  droplet.  Instead,  the  metal  atoms  appear  to  form 
a  weakly-bonnd  “atomic  foam”  in  which  individnal  dopant  atoms  remain  well  separated 
from  one  another  [1,  3]. 

Compntational  stndies  of  metal-doped  He  droplets,  based  on  helinm  density  fnnctional 
theory  (He-DFT)  treatments  of  the  droplet,  snggest  that  the  hrst  solvation  shell  of  He 
atoms  snrronnding  a  metal  atom  dopant  is  likely  to  have  a  density  snbstantially  higher  than 
that  of  bnlk  liqnid  helinm  [4,  5].  For  Mg  dopant  atoms,  as  an  example,  the  atomic  density 
of  the  hrst  solvation  shell  is  predicted  by  He-DFT  stndies  to  be  as  high  as  35  atoms/nm^, 
or  abont  60%  higher  than  the  density  of  bnlk  liqnid  helinm  [4] .  (This  density  is  eqnivalent 
to  the  density  of  solid  helinm  at  a  pressnre  of  abont  130  atm  [6].)  Interactions  between 
the  dense  solvation  shells  that  snrronnd  the  individnal  dopant  atoms  are  thns  thonght  to 
play  a  key  role  in  stabilizing  the  hypothesized  “atomic  foam”  against  recombination  of  the 
metal  dopants. 

A  detailed  investigation  of  the  phenomena  that  contribnte  to  the  stability  (or  metasta¬ 
bility)  of  these  foams  conld  not  only  help  researchers  nnderstand  the  foams’  nniqne  prop¬ 
erties,  bnt  conld  also  shed  light  on  the  physics  that  governs  the  initial  stages  of  metal 
aggregate  formation  in  He  nanodroplets.  This  conld  in  tnrn  help  experimental  scientists 
design  protocols  for  nsing  He  droplets  to  synthesize  metal  atom  aggregates  with  well- 
dehned  compositions  and  strnctnres.  Fonr  recent  reports  [7-10]  of  metal  clnster  formation 
in  He  nanodroplets  demonstrate  the  versatility  and  power  of  nanodroplet-assisted  clnster 
synthesis. 

We  have  developed  a  parallel  code  snite  (powered  by  the  Message  Passing  Interface 
parallel  compnting  protocol)  for  the  stndy  of  the  strnctnre  of  neat  and  doped  He  droplets 
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using  the  variational  path  integral  (VPI)  quantum  Monte  Carlo  (QMC)  method  [12,  13]. 
The  code  suite,  QDROP,  permits  fully  atomistic  simulations  of  these  systems  that  will  com¬ 
plement  the  continuum  He-DFT  studies  cited  above.  Specihcally,  the  VPI  QMC  approach 
mitigates  three  shortcomings  of  the  He-DFT  studies: 

(1)  The  He-DFT  simulations  are  based  on  a  helium  density  functional  that  is  cali¬ 
brated  to  the  properties  of  bulk  superfluid  liquid  helium.  As  we  noted  earlier,  however, 
the  solvation  shells  around  metal  atoms  in  He  nanodroplets  can  have  densities  that  are 
characteristic  of  solid  helium  under  moderate  pressure,  and  it  is  uncertain  whether  the 
helium  density  functionals  typically  used  in  He-DFT  studies  of  metal-doped  He  droplets 
are  reliable  at  these  higher  helium  densities. 

(2)  The  He-DFT  studies  treat  the  metal  atoms  as  inhnitely  massive,  immobile  im¬ 
purities.  A  fully  atomistic  approach  that  treats  the  helium  atoms  and  metal  dopants  on 
an  equal  footing  can  help  us  understand  whether  dopant  atom  zero-point  motions  play  an 
important  role  in  the  formation  of  the  hypothesized  atomic  foams.  Earlier  theoretical  stud¬ 
ies  suggest  that  dopants  solvated  inside  He  droplets  undergo  large-amplitude  zero-point 
excursions  away  from  the  droplet  center  [11],  and  an  accurate  treatment  of  these  large- 
amplitude  excursions  could  be  important  in  a  quantitative  investigation  of  weakly-bound 
dopant  aggregates  inside  He  droplets.  By  artihcially  increasing  the  dopant  atom  mass  in 
our  simulations,  we  can  computationally  suppress  dopant  zero-point  motion;  this  will  allow 
us  to  investigate  the  role  that  dopant  zero-point  motion  plays  in  the  solvation  process. 

(3)  A  fully-atomistic  treatment  allows  us  to  incorporate  three-body  interactions  into 
the  underlying  potential  energy  surface  that  describes  the  doped  He  droplet.  At  the 
elevated  helium  densities  that  characterize  dopant  solvation  shells  in  He  droplets,  it  may  be 
important  to  include  three-body  forces  if  we  are  to  treat  the  solvation  shells  quantitatively. 
Because  the  helium  density  functionals  typically  used  to  study  dopant  solvation  in  He 
droplets  are  calibrated  to  the  properties  of  bulk  superfluid  liquid  helium,  any  three-body 
interactions  that  are  implicitly  included  in  the  parameterization  of  the  functional  represent 
an  “average”  three-body  interaction  characteristic  of  the  liquid.  Three-body  interactions 
could  be  more  signihcant  at  the  higher  densities  found  in  dopant  solvation  shells. 

VPI  QMC  simulations  employ  imaginary-time  propagation  techniques  to  obtain  in¬ 
formation  about  the  quantum  mechanical  ground  state  of  a  many-body  system.  Oper¬ 
ationally,  the  simulations  map  a  quantum  mechanical  V-atom  system  onto  a  classical 
polymer  of  P  “beads”;  each  bead  represents  a  conhguration  of  the  underlying  V-atom 
system  and  moves  in  a  3V-dimensional  coordinate  space.  Neighboring  beads  along  the 
polymer  chain  are  coupled  by  harmonic  oscillator  “bonds”  that  represent  the  portion  of 
the  imaginary-time  propagator  that  is  associated  with  the  quantum  mechanical  kinetic 
energy  operator.  As  P  ^  oo,  the  polymer’s  two  end  beads  attain  asymptotic  probability 
distribution  functions  that  are  equivalent  to  the  V-atom  system’s  ground  state  wavefunc- 
tion,  while  the  polymer’s  interior  beads  attain  asymptotic  probability  distributions  that 
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are  equivalent  to  the  A^-atom  system’s  ground  state  probability  distribution.  VPI  QMC 
simulations  can  therefore  simultaneously  provide  both  the  ground  state  energy  of  the  N- 
atom  system  and  information  about  the  expectation  values  of  coordinate-space  observables 
(such  as  pair  correlation  functions  or  one-body  densities).  VPI  QMC  simulations  can  read¬ 
ily  be  carried  out  on  parallel  computing  platforms,  and  exhibit  nearly  linear  speedup  with 
the  number  of  processors  devoted  to  the  simulation  [14-17]. 

In  Fig.  1  we  show  representative  density  prohles,  computed  using  our  QDROP  code, 
for  Ne-doped  He  droplets  with  V  =  75  to  600  atoms.  This  hgure  shows  the  evolution,  with 
increasing  N,  of  the  solvation  shells  that  surround  the  Ne  dopant.  The  hrst  solvation  shell 
is  already  present  even  at  V  =  75,  but  does  not  reach  its  asymptotic  density  of  roughly 
2.3  times  the  density  of  bulk  liquid  He  until  the  droplet  has  N  =  600  atoms.  The  second 
and  third  solvation  shells  can  be  easily  recognized  in  the  larger  droplets;  a  fourth  solvation 
shell  is  just  beginning  to  emerge  at  N  =  600. 

At  present,  QDROP’s  execution  time  scales  quadratically  with  the  number  of  atoms 
in  the  He  droplet.  A  ten-fold  increase  in  the  number  of  atoms  in  the  droplet  thus  requires 
a  100-fold  increase  in  the  computational  resources  devoted  to  the  problem.  By  using  linked 
cell  lists  [20]  to  exclude  distant  atom-atom  pairs  from  the  computation  of  the  total  potential 
energy  of  the  droplet,  we  can  improve  QDROP’s  scaling  to  be  nearly  linear  with  the  number 
of  atoms  in  the  droplet.  This  will  allow  us  to  carry  out  VPI  QMC  studies  of  droplets  with 
N  ^  5000  atoms  using  roughly  O(IO^)  nodes  of  a  massively  parallel  computing  system  and 
about  one  week  of  wall  clock  time.  Droplets  of  this  size  are  comparable  to  those  employed 
in  experimental  pick-up  experiments  [21]. 

2.2  Molecular  dynamics  simulations  of  droplet  soft  landings 

Classical  molecular  dynamics  simulations  are  a  mature  computational  workhorse  for 
studying  real-time  chemical  dynamics  in  large  systems  (with  more  than  a  dozen  or  so 
atoms).  The  molecular  dynamics  approach  has  recently  been  applied  to  study  the  des¬ 
olvation  of  species  embedded  in  nanodroplets  by  colliding  the  droplet  with  a  solid  sur¬ 
face  [22].  In  these  simulations,  droplets  containing  as  many  as  6000  particles  (typically 
water  molecules  or  argon  atoms)  are  studied  for  periods  of  time  as  long  as  300  picosec¬ 
onds.  This  demonstrates  that  a  molecular  dynamics  study  of  the  soft  landing  of  doped  He 
nanodroplets  containing  as  many  as  5000  atoms  could  be  technically  feasible. 

The  main  drawback  to  using  classical  molecular  dynamics  methods  to  simulate  real¬ 
time  dynamics  of  He  nanodroplets  is  that  these  methods  ignore  both  the  substantial  quan¬ 
tum  mechanical  delocalization  of  the  individual  helium  atoms  in  the  droplet  and  the  quan¬ 
tum  statistical  effects  that  accompany  the  interchange  of  identical  particles  in  a  many- 
body  quantum  system.  While  including  quantum  statistics  in  a  real-time  simulation  of 
5000  identical  atoms  is  nearly  intractable  at  the  moment,  there  are  a  variety  of  strategies 
for  including  the  delocalized  nature  of  individual  atoms  in  a  real-time  classical  dynamical 
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simuulation  of  moderately  quantum  systems  like  liquid  H2  or  liquid  helium.  One  of  these  is 
zero-point  averaged  dynamics  (ZPAD),  which  incorporates  the  effects  of  atomic  delocaliza¬ 
tion  by  treating  individual  atoms  as  smeared-out  probability  distributions  that  move  along 
classical  trajectories.  ZPAD  simulations  have  recently  been  used,  for  example,  to  carry 
out  30-ps  molecular  dynamics  simulations  of  doped  He^  droplets  with  N  ^  100  [22-24]. 

ZPAD  simulations  treat  individual  He  atoms  in  a  droplet  as  frozen,  spherically- 
symmetric  probability  density  distributions  |'?/^(r)p,  obtained  as  the  square  of  a  single¬ 
particle  He  atom  wavefunction  'ip{r).  The  single-particle  wavefunction  is  obtained  through 
an  iterative  procedure  that  involves  solving  the  three-dimensional  Schrodinger  equation 
for  an  individual  He  atom  moving  under  the  influence  of  its  neighbors.  The  ZPAD  method 
traces  its  genesis  back  to  Einstein  models  [25-28]  for  the  zero-point  motions  of  atoms  and 
molecules  in  low  temperature,  highly  quantum  Boltzmann  solids. 

One  computational  advantage  of  the  ZPAD  approach  is  that,  because  the  single-atom 
probability  density  distributions  |'0(r)p  are  computed  “once  and  for  all”  at  the  beginning 
of  the  simulation,  the  underlying  interatomic  interactions  that  govern  the  droplets’  impact 
dynamics  can  be  pre-averaged  over  these  probability  density  distributions  to  give  effec¬ 
tive  interatomic  potential  energy  curves  that  have  zero-point  averaging  effects  built  into 
them  [25].  This  makes  ZPAD  simulations  no  more  expensive  than  conventional  classical 
molecular  dynamics  simulations. 

However,  one  drawback  of  the  ZPAD  approach  is  that  because  the  He  atoms’  proba¬ 
bility  density  distributions  are  dehned  at  the  beginning  of  a  simulation,  these  distributions 
cannot  respond  to  changes  in  the  atoms’  local  environment  that  occur  during  a  simulation. 
As  a  He  droplet  makes  a  soft  landing  on  a  solid  surface,  we  might  expect  the  probability  dis¬ 
tributions  of  He  atoms  near  the  surface  to  temporarily  acquire  some  degree  of  anisotropic 
character,  as  shown  in  Fig.  2.  We  are  developing  a  classical  dynamics  method  that  extends 
the  ZPAD  method  to  incorporate  these  effects  in  simulations  of  He  nanodroplets;  here,  we 
summarize  our  hrst  steps  in  this  project. 

In  the  approach  we  propose,  we  will  allow  the  He  atoms’  single-particle  wavefunctions 
to  vary  during  the  simulation.  Let  t)  represent  the  single-particle  wavefunction  for 

He  atom  number  k  at  time  t.  We  will  expand  this  function  in  terms  of  a  basis  set  of  s 
localized  functions  {^^(u)  :  p  =  l,2,...,s}  centered  on  the  He  atom’s  current  average 
position  {r^): 

S 

t)=Yl  “  ^^k))  ■  (1) 

p=i 

We  then  recompute,  at  each  step  along  the  molecular  dynamics  trajectory,  the  expansion 
coefficients  {c^  :  p  =  1, 2, . . . ,  s}  for  atom  k.  This  is  akin  to  ab  initio  molecular  dynamics 

(AIMD)  approaches  that  carry  out  “on  the  fly”  molecular  dynamics  simulations  by  re¬ 
computing  electronic  wavefunctions  at  each  step  along  the  molecular  dynamics  trajectory; 
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here,  the  He  atom  single-particle  wavefunctions  are  the  analogues  of  the  electronic  degrees 
of  freedom  in  AIMD  simulations. 

By  requiring  the  He  atoms’  single-particle  wavefunctions  to  be  constructed  as  linear 
combinations  of  the  basis  set  functions  {(/)p(u)  :  p  =  1,  2, . . . ,  s},  we  can  provide  the  atoms 
with  flexibility  to  change  their  wavefunctions  as  the  molecular  dynamics  trajecory  evolves, 
but  still  retain  most  of  the  efficiency  of  the  conventional  ZPAD  approach  in  evaluating  zero- 
point  averaged  He-He  interactions.  If  V (r)  is  the  bare  He-He  interatomic  interaction,  then 
the  zero-point  averaged  interaction  for  He  atoms  j  and  k  is  given  by 


(2) 


where  the  angle  brackets  represent  averaging  over  the  two  atoms’  zero-point  probability 
distribution  functions.  Inserting  the  dehnition  of  from  above,  and  letting  = 

we  obtain 


s  s  s  s 


p=l  p'  =  l  q=l  q'  =  l 


(3) 


where  R  =  —  u^-  -1-  S  and  S  =  (r^)  —  {r-)  is  the  vector  that  connects  the  current  mean 

positions  of  atoms  j  and  k.  The  angle  brackets  in  Eq.  (3)  represent  averaging  over  the  two 
three-dimensional  vectors  and  u^i.;  the  quantity  in  angle  brackets  is  therefore  a  function 
only  of  the  vector  S. 

Because  the  basis  set  functions  (/>p(u)  have  functional  forms  that  do  not  change  during 
the  simulation,  the  averaging  process  represented  in  Eq.  (3)  can  be  carried  out  at  a  set  of 
discrete  S  vectors,  arranged  on  a  large,  dense  three-dimensional  grid,  at  the  start  of  the 
simulation.  On-the-fly  computation  of  the  interaction  between  He  atoms  j  and  k  can  then 
be  done  by  interpolating  these  pre-averaged  functions  on  the  S  grid,  and  combining  the 
pre-averaged  functions  with  the  weighting  coefhcients  as  indicated  by  the  summation  in 
Eq.  (3). 

Our  initial  choice  for  the  basis  set  functions  0p(u)  will  be  a  set  of  low-lying  three- 
dimensional  harmonic  oscillator  functions.  Our  approach  for  choosing  the  force  constant 
for  the  oscillator  is  described  in  more  detail  below.  At  each  step  of  the  simulation,  the 
coefficients  {Cp  j}  that  dehne  the  zero-point  motion  for  He  atom  j  will  be  adjusted  to 
minimize  the  He  atom’s  zero-point  energy  in  the  held  generated  by  nearby  atoms.  Using 
a  set  of  harmonic  oscillator  functions  thus  provides  the  He  atom  with  hexibility  to  adjust 
its  own  local  wavefunction  as  the  trajectory  evolves. 

We  plan  to  use  our  modihed  ZPAD  approach  to  simulate  soft  landings  of  doped  He 
nanodroplets,  containing  O(IO^)  He  atoms,  in  which  the  dopant  is  either  an  individual 
Mg  atom  or  the  weakly-bound  Mg2  dimer.  We  will  study  soft  landings  at  surface  impact 
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velocities  that  are  chosen  in  consnltation  with  AFRL  personnel  to  mimic  realistic  exper¬ 
imental  conditions.  For  Mg2  dopants,  we  will  look  for  signs  of  Mg2  dissociation  as  the 
doped  droplet  hits  the  snrface.  The  primary  goal  of  these  simnlations  will  be  to  nnder- 
stand  what  droplet  sizes  and  impact  velocities  might  lead  to  landings  that  are  soft  enongh 
to  minimize  energy  transfer  from  the  droplet  to  the  Mg2  dopant,  and  that  therefore  allow 
the  Mg2  dopant  to  arrive  at  the  snrface  intact  and  withont  snbstantial  vibrational  excita¬ 
tion.  For  simplicity,  we  will  ignore  the  atomic  strnctnre  of  the  landing  snrface  and  treat  it 
as  a  smooth,  impenetrable  wall  with  a  weak  van  der  Waals  physisorption  well  for  the  Mg 
or  Mg2  dopant.  We  will  approximate  the  droplet’s  potential  energy  as  the  snm  of  pairwise 
additive  He-He,  He-Mg,  and  Mg-Mg  interactions;  these  are  averaged  over  the  He  atoms’ 
zero-point  motions  as  indicated  in  Eq.  (3)  to  yield  effective  interatomic  interactions  at  each 
step  of  the  molecnlar  dynamics  trajectory. 

We  also  plan  to  nse  onr  modihed  ZPAD  approach  to  simnlate  the  pick-np  of  Mg 
dopants  by  He  nanodroplets.  We  will  hrst  simnlate  the  pick-np  of  a  single  Mg  atom  by 
a  droplet  to  nnderstand  the  dynamical  processes  that  dissipate  (e.g.  throngh  He  atom 
boil-off)  both  the  relative  kinetic  energy  of  the  atom-droplet  pair  and  the  solvation  energy 
released  when  the  atom  dissolves  in  the  droplet.  We  will  then  simnlate  the  pick-np  of  a 
second  Mg  atom  by  a  droplet  already  doped  with  a  single  Mg  atom;  these  simnlations  will 
complement  the  VPI  QMC  stndies  described  above  by  helping  ns  nnderstand  whether  Mg 
atoms  can  effectively  snrmonnt  any  recombination  barrier  that  arises  from  their  dense  He 
solvation  shells. 

The  force  constant  that  dehnes  the  harmonic  oscillator  basis  set  for  a  given  He  atom 
will  be  obtained  throngh  a  ZPAD-like  approach.  To  test  this  approach,  we  have  applied 
it  to  an  ensemble  of  conhgnrations  that  represent  fragments  of  distorted  hcp-like  solid  He 
crystals.  Specihcally,  we  take  a  13-atom  fragment,  consisting  of  a  central  atom  and  its 
12  nearest  neighbors  in  the  hep  crystal  lattice.  We  scale  the  nearest-neighbor  distance  of 
the  13-atom  fragment  so  that  the  fragment  represents  a  conhgnration  drawn  from  a  lattice 
with  a  particnlar  density.  We  then  randomly  jiggle  the  three  Cartesian  coordinates  of  each 
atom  by  adding  a  random  displacement  of  np  to  ±10%  of  the  nearest-neighbor  distance. 

Eqnipped  with  this  distorted  hcp-like  lattice  fragment,  we  then  assnme  that  each  atom 
in  the  fragment  is  represented  by  a  three-dimensional  isotropic  harmonic  oscillator  single¬ 
particle  wavefnnetion.  We  adjnst  the  force  constant  of  these  wavefnnetions  (keeping  the 
force  constant  for  all  13  He  atoms  the  same)  to  minimize  the  total  (kinetic  pins  potential) 
energy  of  the  central  He  atom.  Once  we  have  minimized  this  energy,  we  then  record  the 
mean  sqnared  displacement  {v?)  of  the  central  He  atom;  this  contains  the  same  information 
as  the  force  constant,  bnt  is  physically  more  meaningfnl. 

Fignre  3  shows  the  distribntion  of  {u^)  valnes  for  1000  distorted  hcp-like  cages  at 
three  different  target  densities:  22.6  nm“^,  27.8  nm“^,  and  34.8  nm“^.  These  densities 
are  respectively  3%,  28%,  and  59%  higher  than  the  density  of  bnlk  liqnid  He,  and  have  been 
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chosen  to  span  a  range  of  densities  that  might  be  present  in  local  regions  of  metal-doped 
He  droplets. 

We  see  from  Fig.  3  that  at  each  nominal  density,  the  range  of  force  constants  (or  mean 
sqnared  displacement  valnes)  derived  nsing  this  procednre  spans  a  relatively  narrow  range, 
and  appears  to  be  represented  by  a  ronghly  Ganssian  probability  distribntion.  This  sng- 
gests  that,  for  onr  modihed  ZPAD  approach,  the  force  constants  that  are  associated  with 
the  harmonic  oscillator  basis  set  nsed  to  represent  the  He  atoms’  single-particle  wavefnnc- 
tion  can  be  qnickly  and  accnrately  estimated  from  the  local  He  atom  density  associated 
with  each  atom’s  immediate  environment.  (This  local  density  can  in  tnrn  be  estimated 
from  a  Voronoi  tesselation  of  the  three-dimensional  space  in  which  the  He  atoms  move.) 
Becanse  of  the  relatively  narrow  range  of  force  constants  associated  with  a  given  nominal 
local  He  density,  it  appears  that  onr  modihed  ZPAD  dynamics  will  be  relatively  insensitive 
to  the  exact  valne  of  the  force  constant  chosen  for  each  He  atom  in  the  droplet;  however, 
this  hypothesis  remains  to  be  tested  throngh  additional  nnmerical  experiments. 
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Fig.  1:  Evolution  of  the  solvation  struc¬ 
ture  in  Ne-doped  He^  nanodroplets  as  a 
function  of  N.  From  left  to  right,  the  den¬ 
sity  prohles  shown  are  for  N  =  75,  150, 
300,  and  600  He  atoms.  The  Ne  atom  is 
held  hxed  at  the  center  of  the  droplet. 


Fig.  2:  A  schematic  depiction  of  two  He 
atoms’  trajectories  in  a  molecular  dynam¬ 
ics  simulation.  Both  atoms  carry  single¬ 
particle  wavefunctions  along  their  trajec¬ 
tories,  represented  by  the  “bulls-eye”  con¬ 
tour  plots  at  selected  points  along  the  tra¬ 
jectory.  One  atom  collides  with  a  repul¬ 
sive  wall  (at  the  bottom  of  the  hgure)  and 
then  with  a  second  atom.  Note  the  slight 
deformation  of  the  single-particle  wave- 
functions  during  the  collisions. 


Fig.  3:  Probability  distributions  for  the 
mean  squared  displacement  {u^)  for  He 
atoms  in  distorted  hcp-like  cages.  The  en¬ 
semble  of  distorted  cages  is  generated  fol¬ 
lowing  a  procedure  described  in  the  text. 
The  nominal  densities  of  the  systems  are, 
from  left  to  right:  34.8  nm“^  (blue),  27.8 
nm“^  (red),  and  22.6  nm“^  (green). 
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